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INTRODUCTION 



Since decades the so-called Rosenbluth separation has been utilized to extract the elec- 
^ ■ tromagnetic form factors of the nucleon from the differential cross section for unpolarized 

2 ! electron nucleon scattering. More recently, polarization measurements have become available 

that provide additional information on these form factors. Surprisingly, substantial incon- 
^ sistencies between these two measurements of the same physical quantity seem to emerge. 

^ ■ To start discussing the problem we introduce the relevant Lorentz invariant kinematical 

! variables. First we have 

^: ^ = Sr = -4SJ <^-^' 

■ where q = k — k' = p' —p is the momentum transfer and M the nucleon mass. For space-like 
O I processes, such as the elastic electron nucleon scattering, = —q^ is non-negative. The 
Q \ second Lorentz invariant variable is the photon polarization parameter 

X. '-ottItti "''^'■^ . = i(«.+r).(p+pr (1.2) 

b : 

The leading (tree level) contribution to elastic electron nucleon scattering is shown as 
Feynman diagram in figure [H In this one-photon exchange (or Born) approximation the 
unpolarized elastic electron nucleon scattering cross section is 

■ Mott e(l + r) ^ • ^ 

in the limit of vanishing electron mass, mg 0, which is well justified in the considered 
kinematical regime. In eq. fll.3p Ge.m are the electromagnetic form factors of the nucleon 
to be further specified later. The Rosenbluth separation to extract these form factors from 
data is to display the ratio of the measured cross section and the Mott cross section at a 
given momentum transfer, i.e. r as function of e. After multiplicating this function with 
e(l + r), the intercept and slope yield Gm and Ge, respectively. 

However, since the relative contribution of Ge to the cross section quickly decreases with 
growing momentum transfer, this type of analysis becomes increasingly difficult already 
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FIG. 1: Tree level loop contribution to the electron (momenta k and A;') nucleon momenta p and 
p') scattering with a local nucleon two-photon vertex. 

at ~ IGeV^. Yet the resuhs of this separation have turned out consistent with the 
assumption that the ratio 

m^)-^. (1.4) 



where /Xp = Gm{0) is the proton magnetic moment, approximately equals unity []|. More 
recently, it has become possible to directly extract this ratio from polarization observables 0] , 
thereby avoiding the potentially uncertain separation technique. Most surprisingly, these 
direct measurements exhibited a linear fall 

R{Q^) ^ 1 - 0.13 (g^[GeV]2 - 0.04) (1.5) 

indicating an eventual root at about ~ lOGeV^. Such a behavior was already suggested 
as early as 1973 within semi-empirical fits to existing data. More recently this structure 
has been predicted within a chiral soliton model study j^. To resolve this puzzle, the 
Rosenbluth separation has been repeated with significantly improved precision [51]. For this 
technique to be operative, it is important that the data are indeed consistent with this linear 
relation [6|. Not only is this indeed the case but also the previous result i? 1 is reproduced. 
Hence we face the paradox situation that two distinct methods to experimentally determine 
a fundamental nucleon property yield inconsistent results 0| • 

A possible resolution could be that contributions to the cross section that stem from two- 
photon exchanges but have been omitted so far, are amplified with increasing Q^. In turn, the 
Rosenbluth separation yields modified form factors that significantly deviate from the ones 
that are defined via a one photon exchange jl]. On the other hand, symmetry properties 
require that the two-photon exchange contributions alter the (linear) dependence on the 
photon polarization parameter e [9|. Yet the data are consistent with but not restricted to 
this linear dependence [13] . It is widely believed [ll[ that the two-photon mainly effects the 
Rosenbluth method but are negligibly small for the polarization process. We will therefore 
focus on the former. 

There are two types of two-photon exchange Feynman diagrams that we display in fig- 
ures 12] and [31 The box diagrams in figure [2] are essentially iterations of the one-photon ex- 
change shown in figure [H Any estimate of these box diagrams require additional information 
about the off-shell behavior of the photon-nucleon vertex. Furthermore the intermediate 
baryon is not restricted to be the nucleon, it could be any nucleon resonance that possesses a 
sizable electro-production potential. Many assumptions and modeling of baryon properties 



enter the computation of these diagrams [8|, [12|, [13|, |14|. |15| . For a recent review on both the 



experimental and the theoretical situations see ref. [iil,[ii. 
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FIG. 2: One loop contribution to the electron (momenta k and k') nucleon (momenta p and p') 
scattering within a typical box diagram. In the hadronic picture, the intermediate baryon can 



be any nucleon resonance [15|], while in the quark picture the photon nucleon vertex can, e.g. be 



related to generalized parton distributions 111 ]. 



The second type of diagrams shown is in figure [3] and has the photon coupled to the 
nucleon at a single vertex. Such diagrams are curious because they do not appear in simple 
Dirac theories of the nucleon which are linear in the covariant derivative. Nevertheless, such 
diagrams are not unknown in hadron physics. In particular, the ir^ —>■ 77 decay induces 
such a vertex as shown in figure HI In that case the intrinsic structure of the 7r°77 vertex 
is dictated by the quark triangle diagram, i.e. the axial anomaly. So we may imagine the 
nucleon coupling to an (off-shell) pion via a Yukawa interaction and the pion subsequently 
decaying into an electron-positron pair as in figure HI As we will observe, this process is 
negligible since (after renormalization) this diagram turns out to vanish in the limit rrie — *• 0. 
The reason is that in the interaction Lagrangian the two photons couple to the derivative of 
the pion field and, when computing the Feynam diagram, an integration by parts produces 
a factor mg. Hence the single pion exchange cannot produce a significant contribution to 
electron proton scattering. However, in chiral model multiple pion exchanges cease to have 
derivative couplings to the two photons and are hence not necessarily suppressed when 
nie — >■ 0. 

The main purpose of the present investigation is to point out that this anomaly induced 
process has a considerable affect for two photon contributions in the Rosenbluth analysis 
and that this process has not been considered previously. Of course, this process by itself 
cannot fully explain the discrepancy to the polarization analysis. 

The required anomaly contribution to the local process shown in figure [3] can be perfectly 
studied within chiral soliton models for baryons. In these models the chiral field U not 
only is the non-linear representation of the pion field but also describes the nucleon as 
a (topological) soliton excitation. In these models we may understand the two-photon 
exchanges shown in figure [3] as the coupling of the nucleon's pion cloud to the electron 
through the anomaly. We will particularly compute the diagram in figure [3] within the 
Skyrme model and focus on the anomaly contribution which is unique because it reflects a 
QCD property. Diagrams that have two photons coupled to the nucleon at a single vertex 
are unique to chiral models because they are formulated in meson degrees of freedom. There 
are also non-anomalous local two-photon couplings. These are quite model dependent. 
Therefore we will here only consider the one induced by the non-linear a model as an 
example and observe that its contribution vanishes as rric — * 0. 

The soliton model successfully describes many properties of baryons and baryon res- 
onances. Though quantitatively the agreement with empirical data is on the 0(1/Nc), 
i.e. 30%, qualitative aspects are well explained, as has been recently reviewed in ref. [iTJ. 
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FIG. 4: One loop contribution to the anomalous decay of the neutral pion into a electron positron 
pair. Left panel: loop diagram from eq. (j3.3p . right panel: contribution from the local counterterm 
eq. (j3.4p . There is an analogous contribution from the loop diagram with the external electron legs 
exchanged, cf. figure [3l 

This makes this model a perfect candidate to explore the above posed problem. The use of 
model calculations is the more unavoidable since in due time lattice calculations for (off- 
shell) two-photon form factors will not be available. Furthermore we recall that in soliton 
models two-photon interactions have previously been successfully considered to compute 
static nucleon polarizabilities [l^. Of course, the current problem goes beyond treating 
these interactions in a static framework. 



II. THE MODEL 

We consider the two-flavor Skyrme model for baryons as the simplest of many models 
that support the soliton picture. In terms of the non-linear representation for the isovector 
pion field vf 

(2.1) 

where r is the vector of Pauli matrices, this model is defined by the Lagrangian 

C = ^tr (d.Ud'^U^) + — ^tr ( [d,U, d^^] [d>^U, a^f/t] ) + ^tr (f/ + f/t _ 2) . (2.2) 

The leading contribution is the non-Iinear-cx-term, which is supplemented by the Skyrme 
term which contains four derivatives on the chiral field and is required to eventually stabilize 



U = exp 



.TT ■ r 
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the soliton. Also the chiral symmetry breaking pion mass term is added. The model parame- 
ters that are determined from meson properties are the pion decay constant = 93MeV and 
the pion mass = 138MeV. On the other hand, the Skyrme term coupling esk ~ 4.0 may 
vary within a certain regime imposed by reproducing nucleon static properties reasonably 
well. 

The hedgehog configuration that builds the soliton and carries unit baryon number reads 



Uo{r) = exp [ir ■ rF{r)] 



(2.3) 



where r = |r|. The chiral angle, F{r) is determined from the stationary conditions that result 
from eq. (12. 2p subject to the boundary conditions F{0) = n and F{oo) = 0. To generate 
baryon states with good spin and isospin we subsequently introduce collective coordinates 
A G SU{2) 

f/(r,t) =A(t)f/o(f)At(t) (2.4) 

that parameterize the spin-fiavor orientation of the hedgehog and quantize them canonically. 
This turns the time variation of the collective coordinates into the spin operator. 



J = a^[Uo] tr 



At 



dA 
~dt 



(2.5) 



where the moment of inertia, a^[?7o] is a functional of the classical hedgehog field. Its 
invariance under combined spin and isospin rotations yields the isospin operator as 



-DabJb with Dab = -ti [ra An A'^] 



(2.6) 



In this relation Dab refers to the adjoint representation Dab of the collective rotations. 

The nucleon wavefunction are Wigner-D functions of the collective coordinates in the 
spin J = ^ and isospin I = \ representation 



{A\J = I,t,s) 



2J+1 



nl/2 



Di:\A) 



(2.7) 



±i and t 



Here \s,t) represents a nucleon state with spin and isospin projections s 
±i, respectively. The identification of total spin and isospin originates from the hedgehog 
structure upon which rotations in coordinate and iso- space are identical. 

In practice we obtain operators for observables that are expressed in terms of the col- 
lective coordinates and their time derivatives. We use eqs (12. 5p and (12. 6p to write the 
latter as operators in the space of the collective coordinates and sandwich them between 
the wavefunctions eq. (12.70 . The matrix elements are finally obtained as integral over the 
collective coordinates which are most conveniently evaluated in terms of their Euler angle 
representation. Later we will particularly require 



{s,t\Dab\s\t') 



■{s,t\IaJb\s' ,t') and 



{s,t\D^aD^,\s' ,t') = ^-{s,t\5ab\s\t') . (2.8) 



These results and techniques are well established in soliton models and we refer to reviews, 
^■9- for derivation and further details. These collective coordinate matrix elements 



come together with the matrix elements that emerge from the spatial dependence of the 
soliton. We will return to their computation in section IV. 
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III. ONE AND TWO PHOTON INTERACTIONS 



We obtain the minimal photon Skyrmion interaction by gauging the Lagrangian. Later 
we will also comment on non-minimal interactions. For the local part, eq. ( 12. 2p this is 
straightforwardly accomplished by replacing the partial derivatives with covariant ones: 



d^U D^U = d^U - teA^ 



Q,U . (3.1) 



Here Q = T3/2 + I/6 and are the quark charge matrix and the photon field, respectively. 
Substituting this prescription not only yields the single photon vertex to the nucleon (repre- 
sented by the soliton) but also higher order interactions, in particular two-photon vertices, 
the so-called seagull terms. For example, the non-linear-cr-term in eq. (12. 2p gives 

^(gauged) ^ :^t^(5^f;5A.f;t)_,:|^^^tr(g[f/t9^f/ + f/9^f/t] 



Q,U 



g,f/M. (3.2) 



The situation is slightly more complicated for the non-local Wess-Zumino term [19( that 
we do not make explicit because it does not contribute to pure hadronic objects in the two- 
fiavor Skyrme model. This term encodes the QCD anomaly and yields one and two-photon 
couplings to the chiral field when gauged with respect to the corresponding U{1) group. The 
techniques to compute these couplings are based on a trial and error scenario to obtain a 
gauge invariant quantity. These techniques are widely described in the literature [l^, 2^ 
and the result for the gauged Lagrangian is 

^wT"" = Y^e'^'''''^{A^tr[Q[U^dMU^d,UU^d^U-UdM^Ud,U^Ud^U^]'^ 

+zeA^d,AptT (2Q^ [U^dM - Ud^U^] + Qd^UQU^ - QUQdM^^ } , (3.3) 

wherein we substituted the physical value of three color degrees of freedom. Most interest- 
ingly this contribution to the action generates a vertex for the neutral pion to anomalously 
decay into two photons via the expansion U = 1 + if ■ n/ f-^ + . . . of eq. (12.11) . Taking 
into account the QED coupling to the electrons, eA^'ife'j^'^e, this then describes the decay 
e~^e~ via the Feynman diagram displayed in the left panel of figure |H This process 



was exhaustively discussed in ref. (21| together with its generalization to 77 ^ /i"*"/! , etc.. 



As a matter of fact, this loop diagram is ultra-violet divergent and induces the counterterm 
Cct. = ^x(A)^e7'^75^etr (2g2 [U^d^U - Ud^U^] + Qd^UQU^ - QUQd^U^^ , (3.4) 

where a = e'^/An = 1/137 is the QED fine structure constant. Furthermore x is a divergent 
coefficient^ such that the sum of the two diagrams in figure |4] is finite. Furthermore A 
refers to the normalization scale that enters to properly reproduce the dimensions of loop 



^ In ref. 21[ counterterm coefficients were independently introduced for the terms in cq. (j3.4[) . However, 



gauge invariance enforces them to appear in exactly the displayed manner. 
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integrals in dimensional regularization with J dH — > ^ J d^l. We separate the finite, but 
renormalization scheme dependent part according to 

Xfin(A) = 6 (^^^ - 7 + ln47r^ - x(A) , (3-5) 

where 7 = 0.577. . . is Euler's constant. We consider this renormalization as effective mod- 
eling of an eventual off-shell form factor for the 7r''77 vertex 22|]. After all, classically there 



is no direct interaction as in eq. (13. 4p between electrons and hadrons. On the other hand, a 
vertex form factor that deceases with the photon momentum renders the diagram in figure H] 
finite and thus does not induce such a direct (and unphysical) electron pion interaction. 

The resulting decay width is most conveniently presented in terms of the ratio with 
respect to the decay into two real photons, 

I>!^.^v?iEi|^(^,|., (3,) 

i (vr" '-f'-f) Svr^m^ t, 

2 

where ^ = The complex amplitude A[^) has a complicated representation in terms of 
Feynman parameter integrals [2l[. Here it suffices to remark that the renormalization scale 
dependence emerges only through its real part, 

ReAiO = Xfin(A) - 61n^ - [in {e)T + 2 [3 - 21n(2)] In (e) + 1(0 (3.7) 

where the reminder, A{C,) is independent of A and finite as rrie — > 0. In this massless limit 
the ultra-violet finite imaginary part 

\n^A{0 = -^L= In + Ve^) (3-8) 
also diverges logarithmically. 



The empirical datum (6.3 ± 0.5) x 10 ^ 2J] for the ratio, eq. (13. 6p is reproduced in the 



range —24 < Xfin(A) < —10 when the renormalization scale in the second term on the right 
hand side of eq. (13.71) is set to A = IGeV. We will adopt that range when investigating 
two-photon processes in the nucleon sector. 

A further remark on the limit rrie — is in order. Eventually we want to assume this 
limit when computing the cross section for electron nucleon scattering since it considerable 
simplifies the kinematics and, of course, is physically meaningful because the energy scales 
that are involved in this scattering process are huge compared to the electron mass. Obvi- 
ously the width, eq. (13. 6p vanishes in that limit. This can be easily understood: Wfj_ is a 
total derivative in the one pion approximation. When shuffling this derivative to the elec- 
tron axial current to which the photons in the loop couple, a factor rrie is produced. So this 
decay goes together with a helicity flip of the electron. Hence our renormalization condition 
prevents us from taking the limit rrie —>■ in the loop. As can be observed from eq. (13. 7p 



the loop itself actually produces a double-logarithmic divergence [21 



IV. NUCLEON FORM FACTORS 

To compute the transition matrix elements for elastic electron nucleon scattering we 
consider eqs. (13. 2p and (13.30 as perturbation and couple the photons to the electrons as 
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indicated in figures [T] and [31 This procedure is standard in QED. However, we also need 
to compute the nucleon matrix elements for the rotating hedgehog configuration, eq. (12.41) . 
These matrix elements are commonly parameterized in terms of form factors. If we extract 
the terms linear in and write it as C^^^ = eJ^A^ the corresponding matrix elements 
introduce Dirac and Pauli form factors via 



{N{p')\J,\N{p)) = U{p') 



U{p) , (4.1) 



where and p'^ are the on-shell momenta of the initial and final nucleons and p'^— p^ 
is the momentum transfer. The above definition is the standard Lorentz covariant parame- 
terization of the matrix elements of the conserved electro-magnetic current, in which U (p) is 
the nucleon Dirac spinor. Note that M is just a parameter in this decomposition and refers 
to the actual nucleon mass rather than the model prediction. It is convenient to introduce 
"electric" and "magnetic" (so called Sachs) form factors 

= F,{Q') - ^F,iQ') , GMiQ') = F,{Q') + F,{Q') , (4.2) 

that show up in the differential cross section, eq. (11. 3p . By pure definition the above form 
factors concern one-photon couplings to the nucleon. Unfortunately, they cannot easily be 
accessed from data because nature does not terminate at that order of perturbation theory. 
Similarly the two-photon couphngs that we extract from eqs. (13. 2p and (13.31) . 

6^1 = e^A^A^S and = eh'^'""' A^d^A.W^ (4.3) 

respectively, also define form factors. While there is only a single and simple Lorentz struc- 
ture for S, 

{N{p')\S\N{p)) = U{p')S^,M)U{p) , (4.4) 
the anomaly requires VF^ to be an axial vector with the decomposition 

{N{p')\W,\N{p)) = U{p') [i,Fa{Q') + q,F,{Q') + za.^q'' Fe{Q')] i,U{p) . (4.5) 

There is no simplification or relation between these form factors as for the ordinary axial 
current because no conservation law applies to VF^. Yet we will see that for unpolarized 
scattering, that concerns the Rosenbluth method, only the first form factor. Fa contributes. 

In the soliton model a major task consists in computing the momentum dependent form 
factors in eqs. (14.11) . (14. 4p and (14.50 . In principle an additional collective coordinate that 
parameterizes the position of the soliton must be introduced in eq. (12.41) . Its conjugate mo- 
mentum will be the linear nucleon momentum. This yields quite a simple recipe to handle 



the linear momentum part of the matrix elements [25|: just take the Fourier transformation 
with respect to (minus) the momentum transfer of the coordinate dependent factors in the 
decomposition of the current operators after substituting the soliton configuration. Essen- 
tially we will have to fold the radial functions in the currents by spherical Bessel functions 
associated with angular momentum of the multiply angular structure. In general we may 
choose any frame to do this calculation. However, it turns out that the Breit frame with 

p = —p' = — and g° = (4.6) 
2 
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is particularly suited not only because it properly reflects the zero energy transfer onto an 
infinitely heavy (large Nq) soliton but also because it directly connects the electric form 
factor, Ge and the magnetic form factor, Gm to the time and spatial components of the 
electro-magnetic current, respectively. In this frame the incoming and outgoing nucleons 
evidently have the same energy E = = p'^ = a/ + Q"^ /A. For baryons with spin ^ we 
the find the Sachs form factors from the matrix elements 

{N{p')\j\Q)\N{p)) = 2MGe{Q'){s',\ss) 

{Nip')\.riO)\Nip)) = -2tGMiQ'y^'qHs',\Sk\ss), (4.7) 

where S is the nucleon spin operator. 

So far we have treated the model in a non-relativistic fashion that restricts the energy 
range to be reliably considered below the nucleon mass. Using the techniques and results of 



refs. |26l . 1391 ] we may extend the nucleon form factor calculation to larger momenta by the 



transformation 



Ge{Q') — 7^'""G^ and G^s{Q') T^^^^Gm , (4.8) 

where 7 = yjl + r is the Lorentz boost factor. Essentially this transforms the from factors 
from the non-relativistic to a relativistic frame. It is worthwhile to note that operator 
ordering ambiguities in quantizing the linear momentum are mitigated by choosing the 
Breit frame as starting point. This is so because we have p'^ = p''^ for this special frame. 



The original study [26[ is based on the Lorentz boost and suggests to put n^; = and 
hm = 1- For more insight it is instructive to reflect on the nature of the transformation, 
eq. (14.81) . Most evidently the momentum interval [0,4M^] of the rest frame is mapped onto 
the space-like momenta in the Breit frame. While small momenta are almost unaffected, 
the form factors at infinity in the Breit frame are obtained from those in the rest frame at 
= 4M^. Even though the latter may be small, there is no general reason for them to 
vanish. In particular, this implies that the form factors do not match the empirical dipole 
form unless ue = um > 2. Thus the values ue = um > 2 are also frequently adopted 
because they are strongly motivated by regardin g th e baryon as a cluster of particles whose 
leading Fock component is a three particle state l4l[ |. In any case, the large behavior 



is not a profound model result but merely originates from the boost prescription and thus 
mainly reflects the kinematical situation. We will henceforth assume ue = um = 2 and 
similarly riwz = 2 for the form factors in eq. (14.51) . We note that there is an additional 
ambiguity in the choice of the mass parameter in the Lorentz factor 7. We take M to be 
the nucleon mass, yet from the point of view from an 1/Nc expansion one could equally 
well argue for the soliton mass which is about 50% larger. Again, this does not significantly 
affect qualitative results. 

In figure Owe show the resulting form factors Ge and Gm and compare them with data. 
The data for Gm are obtained according to the Rosenbluth method while the ratio Gm/Ge 
is taken from polarization measurements that are assumed to be more robust against the 
two-photon contamination. As usual we display these data normalized to the dipole form 
factor, 

Gd{Q^) = o • (4.9) 

(l + QV0.71GeV2)' ^ ^ 

This figure clearly demonstrates that soliton models are able to reproduce the gross features 
of the empirical form factors. This is particularly the case for the linear fall of the ratio 
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FIG. 5: Proton (left) and neutron (right) form factors computed in the Skyrme model as a function 
of momentum transfer after applying the boost. These data have been generated with esk = 3.8. 
The magnetic from factors are additionally normalized with respect to the predicted magnetic 
moments = 6*^(0) ^ 2.33 and Ur, = GfA O) ~ —1.99. The dashed line represents the empirical 



fit, eq. ([13]). Data are from refs. [27|, M, M, M, M, M, m,\M,\35\ 
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FIG. 6: The influence of the non-minimal photon coupling, eq. (I4.10p . on the proton form factors 
for two parameters Lg = 0.0069 (left) and Lg = 0.0045 (right). All other model parameters are as 
in figure [5l 



R{Q^), cf. eqs. ([Ill) and ([LS]). Taking the point of view, that two-photon corrections are 
small for the polarization method ll|] it is hence our task to explain the Rosenbluth cross 
section within the model. 

Nevertheless, deviations of the model predictions from the actual data are apparent. 
Model modifications can improve the agreement with data. For example, it is known from 
chiral perturbation theory studies on the pion radius [s^ that non-minimal photon couplings 
as in 

Cg = -zLgA^^tr [Q {d^Ud.U^ + d^U^dM)] , (4.10) 



are mandatory to correctly reproduce the pion electromagnetic properties. In eq. (14.101) A^'^ 
is the photon field strength tensor and Lg «i 0.0069 is adjusted to the pion radius [37|. This 
term may be understood as resembling the contribution from (short distance) vector meson 
fields that have been integrated out when approximating the effective chiral theory by the 
Skyrme model [s^. In figure we show the proton form factors when the non-minimal 
electro-magnetic coupling, eq. (I4.10p . is incorporated in the the electromagnetic current. As 
can be seen from the right panel of that figure, a moderate adjustment of the additional 
parameter indeed leads to excellent agreement with data. This is particularly the case for 
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the linear decay of R. Being a total derivative, eq. fl4.10p does not affect the form factors at 
zero momentum transfer. 

The reader may also consult ref. 3^ for a more thorough investigation in a vector meson 



soliton model that strongly supports the above statement that soliton models provide a fair 
account of the nucleon form factors, even at large momenta. That investigation also shows 
that the just mentioned deviation from the data can be removed by fine tuning the model. 
The Lorentz boost, eq. (14. 8 p is crucial to gain that agreement. As explained in ref. 3^ 



the strong decrease of the ratio R{Q^) then emerges naturally in chiral soliton models as it 
basically stems from the isospin being generated from a rigid rotation in flavor space. In any 
event, we are mainly interested in whether or not the two-photon exchange coupled via the 
anomaly significantly contributes to the Rosenbluth cross section. To answer this question 
qualitatively, no further fine tuning of the model to reproduce the from factors in detail is 
required. 

As will be shown in the next chapter, only the form factor Fa contributes to the inter- 
ference with the one-photon exchange. For its computation in the Breit frame we consider 
the spatial components of in eqs (13.31) and (13. 4p . 

{N{p')mN{p)) = ±^X^ [Ho{Q')a. + H^iQ') (a. - 3q,q ■ a)] x ■ (4.11) 

Here x refers to the two-component nucleon spinor. Doting this matrix element into q as 
well as averaging the q directions yields two relations between and the form factors in 
eq. (14. 5 P from which 

FAiQ') = [HoiQ') + H2{Q')] (4.12) 

is extracted. The two signs refer to proton and nucleon, respectively. In the next step we 



compute the left hand side of eq. ( 14. lip in the soliton model with the techniques of ref. 25 
We substitute the rotating hedgehog configuration, eq. (12. 4p into the expression for Wi that 
we extracted from eqs (13.30 and (13.40 and take matrix elements between nucleon states. 
They elements are straightforwardly evaluated with the help of eq. (12.80 and by noting that 
A^QA = gl + \D^iTi. Finally we encounter the Fourier transforms of the chiral angle in the 
form 



HoiQ^) = MM drr^ 
Jo 

/•oo 

H2iQ^) = MM drr^ 
Jo 



dF 2 ^ ^ 
— — I — smFcosF 
dr r 

dF 1 ^ ' 

— smFcosF 

dr r 



jo{\q\r) 

J2{m, (4.13) 



where je{z) denotes the spherical Bessel functions associated with orbital angular momentum 
i. Once these momentum dependent functions are computed they are subject to the boost 
transformation, eq. (14.80 with nwz = 2. The resulting form factors Hq and H2 are displayed 
in figure [71 

As long as we omit time derivatives of the collective coordinates, the spin operator, 
eq. (12. 5p does not explicitly enter. In any event, it can only occur in Wq which is the only 
component that contains a time derivative. Direct evaluation of (Wq) gives zero for this 
nucleon matrix element after hermitionizing eventual ordering ambiguities. 
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FIG. 7: Nucleoli axial form factors that contribute at the level of the two-photon exchange to the 
elastic electron nucleon scattering. They are computed from eq. (|4.13|) with the boost, eq. (|4.8|) 
applied with nwz = 2 and normalized to the dipole form eq. (|4.9p . Again we adopted esk = 3.8. 



V. RESULTS 

We mainly intend to point out that there is an effect of the anomaly that reveals itself 
as relevant two-photon contributions to electron nucleon scattering processes. Hence we do 
not attempt any fine tuning of parameters. 

To estimate this effect we merely have to substitute the matrix element of Wfj_ in the 
form factor decomposition, eq. (14.51) into the Feynman diagrams of figure [31 We may then 
formally write the corresponding matrix element for the electron nucleon scattering as 



xU{p') [i^Fa{Q') + qMQ^) + za^^q^FEiQ')] 75 U{p} . (5.1) 

Again is the momentum transferred from the electrons (represented by the spinors u{k) 
and u{k')) to the protons via the two photons. The above parameterization is general for 
couplings via axial currents that enter here because of the e-tensor in the Wess-Zumino 
term |9|. Henceforth we will no longer make explicit the dependence on the momentum 
transfer = —Q^, which is negative for this scattering process. 

The effects of the photon-electron loop are contained in the form factors Wi{q^). We 
will discuss the relevant pieces thereof later. We are mainly interested in the (unpolarized) 
interference with the one photon exchange whose transition matrix element is given by the 
electro-magnetic nucleon form factors, c/. eq. (14.11) . 



= i^u{k')Yu{k)U{p) 



U{p) 



2M 



U{p) , (5.2) 



where in the second equation we employed the Gordon decomposition. We sum the interfer- 
ence contributions over electron and proton spins. This computation is significantly simpli- 
fied by using of momentum conservation and the fact that the spinors U{p) and u{k) obey 
free Dirac equations. Only a few structures contribute when we average over polarizations. 
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The intrinsic parity of 75 requires four additional 7 matrices under the Dirac sum. Two of 
which stem from U{p)U{p) = {f)+m)/2 and the equivalent construction from p' . A third one 
comes from the 7^ in the one photon exchange, eq. (15.21) . Hence the forth must originate from 
the from factor decomposition in eq. (15. ip . This directly shows that the contribution from 
the one-photon exchange to the interference is linear in the magnetic form factor Gm- Since 
we also have a Gordon-type decomposition U {p')a^^q'^'y5U (p) = iU{p') {p^ + IbU^p) the 
Dirac structures associated with neither Fp nor Fe satisfy this criterion. A similar argu- 
ment, of course, holds for the electron form factors wi and W2- Finally, the intrinsic parity 
of the e-tensor enforces the sum over polarizations to be antisymmetric under p ^ p' (or, 
equivalently k ^ k'). Up to overall constants, these considerations determine the final result 

M;M^^^ = ^-^wsFaGm [{k ■ pT - {k ■ pf] 

spins ^ 



= 128Tra^ WsFaGmM^^t{1 + t)^-^^ (5.3) 

that formally shows a deviation from the linear e-dependence found in the one-photon ex- 
change approximation, eq. (II. 3p . We have written this equation such as to make explicit the 
dependence on ^(1 — e)/(l + e) as required by general properties and consistency conditions 
for the two-photon interaction j^. 

If Wfj, were a total derivative, Fp would be the only non-zero from factor and thus the 
interference would vanish in the limit rrif, — »• 0. On the other hand, the term of the hadron 
form factor characterizes the —>■ e~e+ decay that we discussed earlier. 

The information about the photon loop is contained in the electron form factor that we 
compute in dimensional regularization. 



ws = -2- I'dy C ' dx 



JO 



1 r-l-y 

6 I dy dx 



x'^m'l — {1 — X — y) yq^ — ie 



4 / X ml — {1 — X — y) yq^ 



In 



A-D V 47rA2 



X(A) 



(5.4) 



We have made explicit the contribution from the counterterm, eq. (13.40 that eventually 
cancels the part that diverges as D — > 4, according to eq. (13.50 . It is interesting to consider 
the leading contribution in the limit of vanishing electron mass, 

^3 ^ ^3 = -21n(^)+7 + ixfin(A)-31n(|!) (5.5) 

and compare it with eq. (13. 7p . First, we notice that the A-dependence is the same, so that 
after fixing the counterterm via the decay vr e~e~^ the model prediction for the unpo- 
larized cross section does not possess any renormalization scale dependence. Second, the 
rrie — singularity is more severe for the decay than for the considered cross section. This 



reflects the fact that the 



In (4 



2 



divergence is buried^ in the electron form factors Wi 



^ In general the denominator in the first integral in eq. 
divergence. However, for the special combination 1 - 
piece drops out. 



[. (|5.4p always yields a double logarithm as the leading 
— a; in numerator of the first integral in eq. (j5.4p this 
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and W2 that do not contribute to the cross section after averaging the polarizations. In turn 
this imphes that for a prescribed amphtude A in eq. (13. 7p . Xfin = 0[lnmg]^. Hence, for 
momenta ^ the counterterm would actually dominate this contribution to the cross 
section if we chose such a renormalization condition^. The discussion of this singularity in 
ref. 22 1 suggests that it arises independently of the high momentum treatment. In ref. 
the authors carefully analyze the (double)-logarithmic singularities in the box diagrams of 
figure O In that case, single logarithmic singularities emerge always while double logarith- 
mic singularities only occur when the momentum transfer roughly equals the mass of the 
exchanged hadronic resoncance. The triangle diagrams (figure [3]) do not contain hadronic 
resonances. Thus the non-existence of double logarithmic singularities in their contribution 
to electron nucleon scattering is expected; even though they may generally emerge as the 
rare 7r° decay exemplifies. 

To facilitate the discussion of our numerical results for the two-photon exchange contri- 
bution from the Wess-Zumino term we identify the tree level cross-section, eq. (11. 3p and 
introduce the reduced cross-section, 



(-) 



(1 



da 



Mott 



(5.6) 



so that 



da 

dn 



Gli + -G\ + ^—w^M^FAGu^r{\^T){\~e^) 

T 71 



ri2 

^ M 



(5.7) 



For small e the predicted correction to the Rosenbluth form is O (e^) and thus small. 
Note, that in the form of the second equation the ambiguities in choosing the powers nx, 
cf. eq. (14.81) cancel within the square bracket. 

In figures [H] and [H] we present our result for (^)^. To mitigate the model deficiencies 
associated with Gm we compare the model prediction 



'^mod(e) 



1 _|_ £l±E 



1 + v^sM^I^ v/TITTT) 



(5.8) 



to the ratio of empirical data 



{da{T,e)/dn) 



R 



{da{T,0)/dn)j 



(5.9) 



By construction, the experimental value {da{T,0)/dQ)R is the magnetic form factor as ob- 
tained from the Rosenbluth method. We recognize from these figures that the additions from 
the Wess-Zumino term to the unpolarized cross section work into the direction required by 
the data. However, they are about a factor five too small for low Q^. For larger it may 



In a fully renormalizablc theory renormalization eonditions are commonly imposed on Green's functions 
with external legs amputated. This would translate to constrain the off-shell amplitude A rather than 
the physical on-shell decay width F via a renormalization condition. 
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Q" = 2.64GeV" 
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FIG. Results for the normalized reduced cross section, eq. 
data 



compared to experimental 

42l |. The full line is our model result, the dashed line is the model result at the one-photon 
exchange level, i.e. Fa = 0. The dotted and dashed-dotted lines are similarly obtained with 
the polarization result for Ge/Gm-, eq. (jl.Sh substituted for the second term of the numerator in 
eq. (|5.8p . The left and right panels distinguish the normalization at e = and e = 1, respectively. 




fall short by an order of magnitude. When we normalize with respect to e = these ad- 
ditions are strongest around the end-point e — 1. One the other hand, this normalization 
point is not very special and we may adopt equally well e = 1, as displayed in the right 
panel of figure [HI In that case the agreement with data occurs to be significantly better, 
yet it is merely a matter of presentation. It also suggests that this two-photon effect would 
be most strongly pronounced around e ~ 0. This is, of course, not the case as can easily 
be recognized by inspecting eq. (15.31) . However, a common procedure in the literature, that 
discusses the two-photon contamination in terms of the quantity 5 defined via 



In) 



5) 



(5.10) 



and that we show in figure [THl suggests otherwise because of the e dependence of the pre- 
factor. Nevertheless such a presentation is interesting as it disentangles the two-photon 
exchange contribution, i.e. the last term in eq. (15.71) normalized to the one-photon contribu- 
tion to the cross section. Experimentally this corresponds to the ratio of the difference and 
the sum of the cross sections for unpolarized electron-proton and positron-proton scatter- 
ing. There are no new data on this separation. However, the existing data 43|, indicate 
that 6 should be negative 1^ . 
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QlGeV-] 



QlGeV-] 



FIG. 10: The two-photon exchange contribution, 6 in eq. (jS.lOp for two values of the photon 
polarization parameter e as a function of the momentum transfer Q^. Also given are the results 
for different finite parts of the counterterm coefficient, eq. (|3.5p . 



Q =4. lOGeV^ 
= 3.20 Gev' 
= 2.64 Gev' 



FIG. 11: The model prediction for the two-photon exchange contribution, 5 as a function of the 
photon polarization parameter e. 



We observe that the anomaly contribution to 6 has a large slope at small momentum 
transfer, while it levels off with increasing Q^. This dependence is uniform as we vary e. 
As a function of e with fixed we find the largest slope of 6 around e ~ 1, as shown in 
figure [TTJ These features are quite different from the contributions of the box diagrams 
(figure [2]) that were estimated in ref. 12| outside of any soliton model. We stress that this 
is not a contradiction, rather the opposite is the case because these results must be added 
in a full computation of the two-photon effects onto the Rosenbluth method. 

After eq. (13. 8p we argued that we set the renormalization scale A = IGeV and considered 
three cases Xfin(A) = —24, —17,-10 that are suggested by the data for the decay 7r° e~e~^. 
We display the corresponding variations for the cross section in figure [101 This uncertainty 
translates into an 10-20% effect for the two-photon contribution to the cross section. These 
minor variations with the uncertainty in the fixing the model parameter from the underlying 
process vr*^ e^e~ is reassuring as it shows that the Inmg effects are not as severe as 
suspected. All results shown in figures E] and [3 refer to the central value Xfin(A) = —17. 

We have also investigated the non-anomalous two-photon vertex that originates from the 
non-linear cr-model, eq. fl4.3p . Since this interaction does not have any derivative operator, 
the corresponding triangle diagram is ultra-violet finite and thus no counterterm is required. 



Yet we find that the corresponding matrix element, M^'^'^ vanishes as 



nip 



0. Hence this 
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interaction gives negligible contribution within the Rosenbluth method, if at all. 

VI. CONCLUSION 

We have performed a model calculation to shed some light on the discrepancies that arise 
from different methods to extract the nucleon electro-magnetic form factors from data. 
These discrepancies are assumed to be resolved by the inclusion of two-photon processes in 
the computation of electron-proton reactions. Here we have focused on the contribution of 
such a process with the least model dependence and fewest assumptions about off-shell form 
factors. This appears to be the anomaly induced two-photon vertex because it actually is 
a QCD property. It naturally emerges from the nucleon pion cloud coupling to the QCD 
anomaly. This particular two-photon exchange contribution to the elastic electron nucleon 
scattering has previously been overlooked presumably because it vanishes in the one-pion 
exchange approximation. However there is no reason for it to be particularly small beyond 
that approximation. For these reasons we focus on this particular process, which of course 
does not fully explain the observed discrepancies by itself. Nevertheless, it is interesting 
to investigate this effect by its own as it has not been considered previously in the context 
of electron nucleon scattering. Of course, there is no reason to assume that this piece 
by itself fully resolves the discrepancy between the Rosenbluth and polarization analyses. 
At face value the corresponding Feynman diagram is ultra-violet divergent and requires 
renormalization. We impose a renormalization condition that reproduces the empirical decay 
width for the process vr" — >■ 77. This is an ad hoc attempt to deal with the (unknown) off- 
shell behavior of the anomalous 7r° — >■ 77 interaction that has been successfully utilized for 
the description of the pion decay. The Skyrme soliton model is a perfect and the simplest 
tool to study this anomalous contribution to electron proton scattering because it provides 
both, the pion cloud picture of the nucleon and a unique description of the pion anomaly 
coupling via the Wess-Zumino term. We do not exclude that more sophisticated models 
might provide more reliable estimates of this effect. 

The so-computed anomaly contribution to the unpolarized cross section has a minor 
effect on the cross section, of the order of a few per cent. This is to be anticipated for 
an order a = 1/137 contamination. Even though this contribution corrects the leading 
order result into the proper direction these corrections are not sufficient to fully explain 
the observed discrepancy. In this context we stress that this anomaly contribution must be 
considered in addition to contributions from the box diagrams in figure [21 Unfortunately, 
their computation is quite model dependent thereby leading to quite some uncertainties, 
in particular at large momentum transfers. Eventually they can be reduced somewhat by 
phenomenological input for t he g eneralized parton distributions from the amplitude of deeply 
virtual Compton scattering |45l |. The studies of ref. [l2| indicate that these box diagrams 
are the most significant for 6 (the two-photon piece in the unpolarized cross section) at 
small e. The anomaly contribution that we have computed enhances S at moderate e so that 
we expect a negative, say about 5%, effect at small and moderate e while at the boundary 
e — > 1 the contributions of both, the anomaly and the box diagrams are compatible with 
zero. 

In general the model contains additional triangle diagram (non-anomalous) type two- 
photon processes like those shown in figure [3] where the two photons couple simultaneously 
to the pion cloud of the nucleon. In the present model it is natural to assume that the 
dominant such process stems from the non-linear a model. We have seen that it vanishes 
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as the momentum transfer is large compared to the electron mass. 

In the next step we will have to investigate the anomaly contribution in the framework 
of the polarization method. In particular the effects of the form factors Fe and Fp, that do 
not show up in the Rosenbluth method, will be of future interest. 
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